The origin of electromagnon excitations in multiferroic iiMnOa 
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Electromagnon excitations in multiferroic orthorhombic iJMnOs are shown to result from the 
Heisenberg coupling between spins despite the fact that the static polarization arises from the much 
weaker Dzyaloshinskii-Moriya (DM) exchange interaction. We present a model incorporating the 
structural characteristics of this family of manganites that is confirmed by far infrared transmission 
data as a function of temperature and magnetic field and inelastic neutron scattering results. A 
deep connection is found between the magnetoelectric dynamics of the spiral phase and the static 
magnetoelectric coupling in the coUinear E-phase of this family of manganites. 
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The coupling between the magnetic and ferroelectric 
order in a diverse set of materials termed multiferroics 
is currently a topic of intense study P, The inter- 
play between these two orders is particularly striking 
in materials where ferroclcctricity appears as a conse- 
quence of spontaneously breaking of inversion symmetry 
of the magnetic ordering. In many such magnetic ferro- 
electrics the spins order in an incommensurate cycloidal 
spiral state Q . The microscopic origin of ferroelectricity 
for this case has been discussed by a number of authors 
[Ij S IIIj leading to a consensus, gcnerically termed the 
spiral mechanism. It relies on the lowering of the en- 
ergy of the antisymmetric Dzyaloshinskii-Moriya (DM) 
exchange in the spiral state by a polar lattice distortion, 
which induces an electric polarization, P oc Q x R, where 
R oc Si X Si+i is the spin rotation axis, and Q is the wave 
vector of the spiral. These ideas have been of central 
importance in the recent discovery of new multiferroic 
compounds. 

Another apparent consequence of multiferroicity is the 
existence of novel coupled magnon-phonon excitations 
called electromagnons 0, • A magnon that gives rise to 
oscillations of electric polarization can be excited by elec- 
tric fields, thereby coupling much more strongly to light 
than the usual magnetic dipole excitation of magnons 
corresponding to antiferromagnetic resonance (AFMR). 
The resulting electric dipole spectral weight has been 
transferred from the phonons down to the magnon fre- 
quency. The dynamic magnetoelectric effects resulting 
from the coupling between spin and polarization waves 
were discussed theoretically at an early stage of the re- 
search on multiferroic materials More recently, Kat- 
sura, Balatsky and Nagaosa (KEN) noted that the mag- 
netoelectric coupling of the spiral mechanism, also gives 
rise to an electromagnon [10|. When the spiral plane 
rotates around Q, so does the induced electric polar- 
ization, which couples this magnetic excitation to elec- 
tric field e of a light wave normal to the spiral plane: 
e||R. The first observation of the electromagnon peak 



for e[|a in TbMnOa with the 6c-plane spiral spin ordering 
(QII&) seemed to confirm this selection rule |7|. However, 
recent measurements on other spiral multiferroics from 
the same family of materials, Euo.75Yo.25Mn03 [llj and 
DyMnOs , showed that this selection rule is violated 
and that instead, the selection rule is tied to the lattice 
(ist . Therefore, the critical questions are the origin of 
the observed excitations, and why the KBN mechanism 
does not seem to apply to these multiferroics. 

The violation of the KBN prediction implies that a 
different dynamical magnetoelectric coupling is respon- 
sible for the appearance of electromagnons in the spiral 
state. In this Letter we report a magneto-infrared trans- 
mission study of TbMnOa and a theoretical model based 
on symmetric Heisenberg exchange that clarify the origin 
of electromagnon excitations in cycloidal i^MnOa. Addi- 
tionally it is shown that this coupling is also responsible 
for the ferroclcctricity in the coUinear magnetic E-phase 
of rare earth manganites , and based on the dynamic 
response of the cycloidal manganites an estimation of the 
polarization in the E-phase is given. 

Single crystals of TbMnOs were grown as described 
elsewhere [ISj . The crystals were oriented by Lane x-ray 
diffraction, and surfaces polished perpendicularly to the 
three crystallographic axes were prepared. Two samples 
were used in these measurements: (1) ab sample for zero 
field experiments of size 3x3 mm^ and 50 /im thick 
along a,b and c, respectively [Phnm setting), and an (2) 
ac sample for magnetic field measurements (H||&) of size 
2x2x1 mm^ for e| |c and then polished to 75 jim for e\\a. 
Transmission measurements were done using a Fourier 
transform spectrometer connected to a magnet cryostat 
with optical access windows in the Faraday configuration. 
Sample (2) allows access to the phase where ferroelectric 
polarization P is changed from the c to the a axis with 
a field of 5 T at 7 K, and the spiral spin configuration 
rotates from the be to ab plane. The transmission spectra 
were then converted to the absorption coefficient oi{v) 
using T{y) = (1 — _R(i'))^e~"^'^^', where I is the sample 
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FIG. 1: (Color Online) (a) Zero field absorption spectra in 
TbMnOa. Main features of spectra correspond to 2 electro- 
magnon peaks at 25 and 60 cm~^ that activate below Tpe, 
and the lowest infrared active phonon just below 120 cm~^. 
Additionally, a broad absorption below the phonon persists 
well above Tjv. (b) Magnetic field dependence (H || b) of 
the low temperature spectra. The frequencies of the electro- 
magnons shift in the high field phase, and their intensities 
increase, (c) Absorption spectrum for the expected polar- 
ization of the KBN electromagnon. The new absorption line 
around 21 cm~^ in the high field phase is magnetic dipole ac- 
tive given that the expected contribution to e (shown in (d)) 
is larger than found in the static measurements [l7| . 



thickness jla ]. 

Figure IDJa) shows the spectra in the far infrared taken 
at zero external magnetic field in TbMnOa. AU the 
features observed have been identified in the spectra of 
Euo.75Yo.25Mn03 [m and are described in the caption. 
We then performed measurements of the far infrared 
transmission in a magnetic field (Hl|6) to test the be- 
havior of electromagnons on different spiral states. The 
spectra taken at 7 K are shown in figure [Ijb). Above ap- 
proximately 5 T the system undergoes a spin flop tran- 
sition from a spiral in the be plane to an ab plane spiral 
yet we observe that the electromagnons remain in 
this configuration e\\a with only a slight shift of their 
frequencies and an increase of their spectral weight. 

It is particularly striking that the absorptions activated 
in the cycloidal phase are only active for light polarization 
e\\a in all systems studied, regardless of the spin plane, 
or static polarization direction or value. This observa- 
tion directly contradicts the prediction of KBN . In- 
deed, the inadequacy of this mechanism is also suggested 
by the large Born effective charge required in the KBN 
model to account for the experimentally observed oscilla- 
tor strength. The observed selection rule clearly indicates 
that account must be taken of the specific characteristics 
of these manganites, such as the crystal structure and 



anisotropic magnetic interactions, in order to find an ac- 
curate description of the electromagnon excitations. 

The KBN prediction for the selection rule of electro- 
magnons when the spiral plane rotates from be to ab was 
tested with measurements in the H||6, e||c, h\\a configura- 
tion. The absorption spectra are shown in Fig. [TJc). The 
spectrum in the high field phase shows clearly an absorp- 
tion feature around 21 cm~^. By fitting the transmission 
spectrum with a model electric dipole active excitation 
we estimate the putative contribution to the static di- 
electric function e as Ae « 1.2 as shown in figure [T[d). 
This change in e is not observed in the static data [17| . 
showing that this feature is not an electric dipole excita- 
tion. Therefore this excitation is a magnetic dipole (/i||a) 
AFMR as expected from the magnetic order. For com- 
parison we note that in the zero field spectra in the e||6, 
h\\c configuration an AFMR around 22 cm~^ is observed 
(not shown). This means that the h\\c magnon rotates to- 
gether with the spin plane thereby changing its selection 
rule to /ijla, consistent with the change of polarization di- 
rection [1| and recent neutron scattering measurements 
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The essential features of magnetic excitations in the 
spiral state can be understood within a model including 
only supercxchangc interactions between Mn spins. 



(1) 



When magnetic anisotropics are neglected, the 4 Mn ions 
in the unit cell of iJMnOa, located at Ri = (i 0,0), 
R2 = (0,i,0), R3 = (iO,i), andR4 = (0,i,|), be- 
come magnetically equivalent (see Fig. [2]). The com- 
petition between the nearest-neighbor ferromagnetic ex- 
change, J = Ji^i+{b±a)/2 < 0; and the antiferromag- 
netic next-nearest-neighbor exchange along the b axis, 
Jb = Ji,i+b > 0, favors the circular spin spiral with 
the wave vector Q\\b, where cos ^ = ^j-, provided that 
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exchange along the c axis, Jc = 
to a 'double spiral' structure with antiparallel spins in 
neighboring ab layers: 



The antifcrromagnctic 
Ji,i+c/2 > 0, gives rise 



Si) = ± I ccosQ-Ti — bsinQ-Ti 



(2) 



where the upper (lower) sign corresponds the ab layers 
with integer (half- integer) z/c. The be plane is favored 
e.g. by the single-ion anisotropy, ^ J2i i^f)'^ with K > 
0, which does not spoil the equivalence of Mn sites (effects 
of other magnetic anisotropics are discussed below). 

The magnon spectrum of the 'double spiral' state 
has one acoustical and one optical branch plotted in 
Fig. Elja) with, respectively, blue dashed and red soHd 
line, as a function of the wave vector k in the co-rotating 
spin frame, in which (S)||z on all sites. For acoustical 
magnons the order parameter L = Si -I- S2 — S3 — S4 
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FIG. 2: (Color Online) One ab layer of the Pbnm unit cell of 
iJMnOa consisting of 4 Mn ions (blue circles) and 4 oxygen 
ions (pink circles) . The displacements of the oxygen ions from 
the midpoints of the straight lines (dashed lines) connecting 
neighboring Mn ions are STi-^-b,2+a ~ (<5x; Sy; Sz), 5ri,2 = 
{-Sx- -3y- —5z), (5ri,2+a = (— (5x; Sy; Sz), and Sri+b,2 = ((5x; 
—5y; —Sz). The labels of the Mn ions in the next layer (with 
z = c/2) are given in parentheses. The displacements of the 
oxygen ions in the next layer have opposite Sz. 



oscillates in the spiral plane, while optical magnons cor- 
respond to out-of-planc oscillations of L. The KBN elec- 
tromagnon is the symmetric superposition of the optical 
magnons with the wave vectors ±Q, while the antisym- 
metric superposition is the AFMR discussed above. 

Next we discuss a mechanism of magnetoelectric cou- 
pling that only involves the isotropic Heisenberg ex- 
change between non-collincar spins and, therefore, is in- 
sensitive to the orientation of the spiral plane. Due to the 
GdFeOa distortion of i^MnOa compounds, oxygen ions 
mediating the supcrcxchangc between nearest-neighbor 
spins in the ab layers are displaced from the straight 
lines connecting two neighboring Mn ions in ab layers 
as shown in Fig. O When an applied electric field shifts 
all oxygen ions by an equal distance along the a axis, the 
exchange constants Ji^2 and Ji^2+a will be changed by 
A J proportional to the applied field, while the exchange 
constants Ji+6,2 and Ji^i,^2+a will be changed by — AJ. 
In other words, due to the alternating rotations of the 
MnOg octahcdra in i?Mn03, a uniform electric field in 
the a direction gives rise to an alternation of the nearest- 
neighbor exchange along the spiral propagation vector 
Q\\b, J oc Jo + AJcos(ko • r), where ko = (0, ^,0). The 
corresponding coupling of spins to the electric field e\\a 
compatible with the Pbnm symmetry has the form: 



5i?a^ [(Si,, 
j 



Slj + b) ■ (S2j + S2j+a) 



(3) 



where the indices 1-4 label different Mn ions in a unit 
cell while j labels cells. 

This interaction couples electric field through the al- 
ternation of the Heisenberg exchange along the b axis to 
the acoustical magnon with kg at the boundary of the 




150 



FIG. 3: (Color Online) (a) The typical magnon spectrum for 
a circular bc-spiral along k = (0,fci,,0), where ki, is measured 
in units of The empty and filled circles mark the posi- 
tion of, respectively, the KBN electromagnon and the electro- 
magnon excited through the rotationally invariant coupling 
Eq.(|3]). (b) Corresponding dielectric susceptibility for the 
(arbitrarily broadened) single-magnon peak (solid line) and 
bi-magnon spectrum (histogram). 



magnetic Brillouin zone of the model Eq.([T]). This ex- 
citation corresponds to relative rotations of the spins Si 
and S2 in the spiral plane, which occur in anti-phase with 
the rotations of the spins S3 and S4 and result in alter- 
nation of the angle between neighboring spins along the 
spiral propagation vector. 

This mechanism of dynamic magnetoelectric coupling 
works only for non-coUinear spins because the first-order 
variation of the exchange energy oc (S.^ • Sj) with respect 
to spin oscillations SSi{t) = {Si{t) — (S^)) ± (S,) is zero 
for coUinear spins. Although interactions of spins with 
e\\b and e||c, similar to Eq.(l3]), are allowed by symme- 
try, they give rise to alternation of the nearest-neighbor 
exchange constants in the directions orthogonal to the 
spiral wave vector Q and do not lead to single-magnon 
excitations by electric field, as neighboring spins along 
the a and c axes are collinear. This explains why the 
polarization of electromagnon in iJMnOa is parallel to 
the a axis independently of the orientation of the spiral 
plane. 

We note that the coupling Eq. ^ induces a static elec- 
tric polarization PeHo in iJMnOs with the collinear an- 
tiferromagnetic ordering of the E-type IJ] , for which the 
scalar product of neighboring spins alternates along the 
6 axis. This fact allows us to estimate the value of the 
polarization in the E-phase using the optical absorption 
data for the spiral phase (the direct measurement of Pe is 
difficult due to absence of single crystal samples) . The to- 
tal spectral weight of the 60 cm"-'^ peak, S « 7000 cm~^, 
corresponds to Pe ~ luC/cm^, in perfect agreement with 
ab initio calculations [19[. For comparison, the spectral 
weight of the KBN electromagnon, calculated from the 
polarization P ~ 8 x 10~'^/iC/cm^ in the spiral phase 
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of TbMnOa, is only « 10 cm"^. The photoexcitation of 
magnons through the magnetoelectric couphng Eq. ([3]) of 
exchange origin, which is about 2 orders of magnitude 
stronger than the relativistic couphng in spiral mangan- 
ites, is much more effective than the KBN mechanism. 

The coupling Eq.([2]), as well as another invariant, 
Ea [(Si - Si+c) • S3 - (S2 - S2+c) • S4], gives rise to the 
photo-excitation of a pair of magnons with the total wave 
vector ko(the so-called 'charged magnons' (l^l)- The cal- 
culated shape of the bi-magnon continuum for the cou- 
pling Eq.Q is shown in Fig. [31^b). Its maximum is lo- 
cated around the electromagnon frequency and its to- 
tal spectral weight for all reasonable parameters of the 
model is less than 10% of the spectral weight of the single- 
magnon peak. This shows that the peak at 60 cm~^ re- 
sults from the photoexcitation of a single zone boundary 
magnon, which can be identified in the ma gno n spec- 
tra obtained by inelastic neutron scattering [2l|. Even 
though charged magnons can also be excited by e||c and 
e\\b through couplings similar to Eq.Q, strong absorp- 
tions in these polarizations are not observed in the ex- 
periments [11 1 . 

Concerning the nature of the weaker 25 cm~^ peak, 
we note that for an elliptical spiral structure the magnon 
wave vector k is not conserved and a number of different 
mechanisms couple high- and low-frequency magnons. 
Thus, the spiral ellipticity resulting from the magnetic 
anisotropy in the spiral plane gives rise to satellite single- 
magnon peaks and bi-magnon continua at the total wave 
vector ko±2Q, ko±4Q . . ., which could correspond to the 
low frequency peak. Furthermore, there are anisotropic 
spin-spin interactions, which are insensitive to the flop 
transition (e.g. involving products SfSj) and couple the 
electromagnon with k = ko to magnons with k = ±Q 
giving rise to low frequency absorption. Which of these 
mechanisms gives the dominant contribution to the low- 
frequency absorption is not clear, but all of them are 
different from the coupling inducing electric polarization 
in the spiral state. 

In conclusion, we have presented experimental data 
and a theoretical model that strongly support the Heisen- 
berg exchange interaction as the origin of the magneto- 
electric dynamics in the family of multiferroics RM11O3. 
Additionally, the fact that these excitations are observed 
both in the neutron measurements and the IR experi- 
ments reflects their origin as single-particle excitations 
of a hybrid electromagnon nature. The surprising out- 
come of this study is that optical data can be used to 
explore properties of competing ferroelectric states: from 
the measured spectral weight of the electromagnon peak 



in the spectrum of spiral manganites we obtained the 
value of the spontaneous polarization in manganites with 
the E-type coUinear ordering. Our results imply that the 
static and dynamic magnetoelectric coupling are differ- 
ent in general and therefore that electromagnons can be 
observed, in principle, in non-multiferroic materials with 
non-coUinear spin orders. This observation opens a new 
avenue of investigation of the dynamic properties of frus- 
trated magnets. 

This work was supported by the NSF MRSEC pro- 
gram under grant DMR-0520471. MM was supported by 
the Stichting voor Fundamcnteel Onderzoek der Materie 
(FOM). 

Note added: During the final preparation of this 
manuscript we became aware of the work of Takahashi 
et al. [22| where they report very similar results in 
TbMnOa. However, they interpret the excitations as 
coming from an electric dipole continuous band of two- 
magnon absorption, which our work in fact rules out. 
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